QED2+1 in graphene: symmetries of Dirac equation in 2 + 1 dimensions 
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It is well-known that the tight-binding Haniiltonian of graphene describes the low-energy excita- 
tions that appear to be massless chiral Dirac fermions. Thus, in the continuum limit one can analyze 
the crystal properties using the formalism of quantum electrodynamics in 2+1 dimensions (QED2+1) 
which provides the opportunity to verify the high energy physics phenomena in the condensed mat- 
ter system. We study the symmetry properties of 2+1-dimensional Dirac equation, both in the 
noninteracting case and in the case with constant uniform magnetic field included in the model. 
The maximal symmetry group of the massless Dirac equation is considered by putting it in the Jor- 
dan block form and determining the algebra of operators leaving invariant the subspace of solutions. 
It is shown that the resulting symmetry operators expressed in terms of Dirac matrices cannot be 
described exclusively in terms of 7 matrices (and their products) entering the corresponding Dirac 
equation. It is a consequence of the reducibility of the considered representation in contrast to the 
3+1-dimensional case. Symmetry algebra is demonstrated to be a direct sum of two gl(2,C) algebras 
plus an eight-dimensional abelian ideal. Since the matrix structure which determines the rotational 
symmetry has all required properties of the spin algebra, the pseudospin related to the sublattices 
(M. Mecklenburg and B. C. Regan, Phys. Rev. Lett. 106, 116803 (2011)) gains the character of the 
real angular momentum, although the degrees of freedom connected with the electron's spin are not 
included in the model. This seems to be graphene's analogue of the phenomenon called "spin from 
isospin" in high energy physics. 

PACS numbers: 03.65.Pm, 73.22.Pr, ll.30.-j, 11.30.Rd 



I. INTRODUCTION 

Graphenei is a single layer of carbon atoms arranged 
into a planar honeycomb lattice. It has attracted a con- 
siderable attention since its fabrication in 2004^1, due to 
its unique attributes being a consequence of charge con- 
jugation symmetry between positive and negative charge 
carriers, which have the internal degree of freedom analo- 
gous to the chirality defined in the quantum field theory^. 
This similarity between condensed matter physics and 
the quantum electrodynamics (QED) makes graphene 
a test-bed for experimental verification of phenomena 
well-known in the particle physics^. For example, the 
Klein paradox can be observed as the electron propaga- 
tion through potential barriers with the graphene playing 
the role of an effective medium^; the Zitterbewegung ef- 
fect restricts the minimum conductivity to the order of 
conductance quantum e? /h (Ref|6|), while the anomalous 
quantum Hall effeclii can be related with the index the- 
orem. 

Such exciting properties and perspectives are a direct 
consequence of the fact that the low-energy properties 
of electrons in graphene can be described by the model 
based on the continuum limit (a — > 0) of the tight binding 
approximatioi:^*^, which obeys a relation formally identi- 
cal to the 2+1 dimensional Dirac equation with the holes 
and the pseudospin states of the A and B sublattices 
being the counterparts of the positrons and the spin, 
respectively^. The remarkable insight into the origin of 
graphene's uniqueness can be achieved by the studies of 
the symmetries of the model. Recently, much attention 



has been paid to the close relation between spatial sym- 
nretries and the existence of the Dirac fermions which 
can be present even in the chenrisorbed graphene with 
the defects distributed with specific symmetryiS. This 
hypothesis have been also verified by means of density 
functional theory (DFT) calculationsii. Moreover, the 
effects of symmetry breaking have been widely studied, 
for example, due to the mass term^ or valley dependent 
vector potentiali^. The latter examples are extremely 
important for applications, since real graphene systems 
always interact with surrounding environnrent which dis- 
turbs its exceptional electronic propertie a^'-'^'^' '— . 

In this context, it is useful to derive the fundamen- 
tal symmetries of the nrodel using the same mathemat- 
ical procedures that are commonly applied to study a 
3+1 Dirac equatioiii^. The main difference between the 
two-dimensional system of graphene's fermions and the 
3+1-dimensional system is that the Dirac equation of 
the latter is based on the four-dimensional irreducible 
representation of Clifford algebra. In contrast, in the 
analysis of the two-dinrensional case one also uses the 
four-dimensional representation, but, as can be easily 
demonstrated, it becomes reducible. The Clifford algebra 
for three-dimensional space-time possesses two inequiv- 
alent irreducible representations which are both two- 
dimensional. Thus, the representation used to study the 
model of graphene is the direct sum of the irreducible 
oires. 

The purpose of this paper is to study the symmetries of 
the 2+1 massless Dirac equation and to clarify the con- 
sequences of the reducibility of the matrices representing 
the Clifford algebra. Understanding the symmetries of 



the underlying dynamics is especially profitable in the 
context of the possibilities of experimental verification 
of 2+1 QED-like phenomena in the condensed matter 
physics as well as tuning of the graphene's properties for 
applications. Due to the significant role of the experi- 
ments concerning the anomalous quantum Hall effect in 
graphene, it is also useful to study the case of 2+1 sys- 
tem in the presence of magnetic field. Moreover, it is the 
only type of interaction that does not disturb the alge- 
braic structure of symmetries of this system and does not 
lead to any qualitative changes in the spectrum (in the 
sense that it keeps the energy gap between particles and 
holes intact). 

In this paper we have determined the algebra of the 
corresponding symmetry operators using the procedures 
elaborated in the case of four-dimensional Dirac equation 
in Refflg. both for free case as well as the case with the 
constant uniform magnetic field. The convenient way to 
treat the latter case is to make use of the supersymmetry 
inherent to the problem. 

The paper is organized as follows: in Sec. II the "on- 
shell" symmetry generators for free massless Dirac equa- 
tion are introduced and the explicit form of its Lie algebra 
generators is found. In Sec. Ill the extended model with 
the magnetic field included is considered and compared 
to the noninteracting case described in Sec. II. The con- 
clusions and perspectives are discussed in Sec. IV. 



II. THE SYMMETRIES OF THE FREE 
MASSLESS DIRAC EQUATION 

We consider the free massless Dirac equation in mo- 
mentum representation: 



7^Pm* = 



(1) 



where 7 matrices are given in the following representa- 
tion: 



/ 
/ 



cr* 
-cr' 



, 7^ = tj^'YYr (2) 



with i — 1,2,3 (although only 7^,7^,7^ enter the con- 
sidered representation of the Clifford algebra). Eq.(IT]) 
implies po = ± \ p \', thus we can write 



7>^==± 1^1(7° ±n,Y) 



(3) 



where nt = pi/ \ p \ with i ~ 1,2, is the unit vector in 
the direction of momentum. We can restrict our analy- 
sis to the case with po =| P \- The minus sign in eq. ^ 
can be accounted for by reversing the sign of the momen- 
tum. Imposing the condition | p 1 7^ one is left with the 
operator: 



L = 7° + n.iY = 



/ 

/ + 71,0-' 



-riia' 




(4) 



Next, we can represent ^ by the Weyl spinors: 

1±75 



* 



*+ 
*_ 



*- 



-^ 



Then eq.([T|) reduces to: 

For any spinor ^ we put (e, e' — ±) 



^t 



- eY 



1 + e'n.a' 



I * 



(5) 



(6) 



(7) 



and ^^ form a new basis given by the following expres- 
sions: 



ei 
63 



2 ' 

2 ^ 



62 = *+ 

64 = *1 



(8) 
(9) 



In this basis the operator L acquires the Jordan block 
form: 



L = 





10 



10 



(10) 



We define the symmetry operators as the matrices S pre- 
serving the null eigenspace of L: 

i* = =^ L5* = (11) 

It is straightforward to compute the general form of S: 



S 



an ai2 

2^11 bu a-12 612 

021 022 

2^21 ^21 2^22 ^22 



(12) 



We can analyze the properties of S in terms of its Lie 
algebra's structure. First, it is worthwhile to note that 
there are two commuting gL(2, C) subalgebras consisting 
of the matrices given in the following form: 



flu ai2 



021 022 





and 





hi 612 



^21 622 



(13) 



Moreover, there is a four-dimensional complex abelian 
algebra A(4, C): 





xu a;i2 



2:21 2:22 



(14) 



Under the adjoint action v4(4, C) provides the represen- 
tation of both gL{2, C) subalgebras. One can conclude 
that the symmetry algebra has the following structure: 



5=(5L(2,C)®gL(2,C))+^(4,C) 



(15) 



where -Q denotes the semi-direct sum. Exphcitely, choos- 
ing the standard basis one finds: 

[Xij , Xmn\ = 

[BijjXmn] = OjmXin (16) 

One can return to the natural basis ^ using the simi- 
larity transformation L = VLV~^ with the matrix: 



gebra in the representation given by eq. © reads: 



V = 



ri- 
0-110 
7i_ 
-10 1 



(17) 



where n± = Tii ± 7712. 

Finally, simple calculations lead to the following result 

for generators of symmetry algebra: 

An = \{l - 7')(1 + 7V) - \n+l\l + 7')7"7- 
A12 = -\l"{l + l'){l 7%') + J"+7°(l + 7')7V 
A21 = -\l"{l - 7')(1 + 7%') - ^ri+7"(l - 7')7V 
Bu = ^(1 + 7')(1 + 7V) - ^^+(1 + 7')7V 
S12 = ^7°(1 - 7')(1 - 7V) + Jri+7"(1 - 7')7V 
B21 = ^7°(1 + 7')(1 + 7V) - J^H7°(1 + 7')7V 
B22 = 1(1- 7')(1 - 7"7') + J"+7"(l - 7')7°7^ 
Xn = ll"{l 7')(1 - 7"7') + Jn+7°(1 " 7')7°7- 
+ 1/(1 _ /)(i + //) + i„_/(i _ Y'h"i^ 

X12 - \l'{l + l'){l + 7%') + Jn_(l + 75)7V 
-1(1 + r-){l 7V) - ^"+(1 + 7')7V 

X21 = 1(1 - 7')(1 - 7V) + \n^il - 7^)7 V 

+ ln+7"(l + 75)7%- + ln_(l - 75)(1 + 7V) 

X22 = J7°(l + 7')(1 + 7"7') + Jn-7"(1 - 7')7V 

-ln+7"(l + 75)7%- + in_7"(l + ^^){l - 7"7-^) 

(18) 

where we denoted 7^ = 7-'^ ± 77^ . 

The explicit form of the generators of the symmetry al- 



An 



i21 



Bll = 



B2I 



Xn 



-''^21 







-71+ 1 



71+ -1 




1 

-71+ 









10 

71+ 0. 

1 -71_ 

-71+ 1 










1 

71+ 



A12 



A22 = 



B12 — 



Bo 



X12 



X. 







-71+ -1 








0" 

71+ 1 


0. 




10- 

-71+ 

0. 


0- 

10 
71+ 0. 




1 71_ 

-71+ -10 




■000 



1 71_ 
.71+ 1 


0- 



0. 



(19) 



Let us note that, due to the fact that the Dirac ma- 
trices 70,71,72 generate a reducible algebra, one cannot 
express all symmetry operators in terms of them onlyii; 
instead, as can be seen from eq. (jlSp . one has to add 73 
matrix to obtain irreducible algebra containing all neces- 
sary matrices. 



III. DIRAC EQUATION IN A CONSTANT 
UNIFORM MAGNETIC FIELD 



We consider the 2+1-dimensional Dirac equation in the 
constant uniform magnetic field B. The general symme- 
try pattern in this case is the result of the following prop- 
erties: (i) the 2+1-dimensional Dirac equation based on 
the irreducible representation of Clifford algebra exhibits 
supersymmetryi^, (ii) the representation of the Clifford 
algebra under consideration is reducible and is the direct 
sum of both inequivalent irreps. 

We include the electromagnetic field through a min- 
imal coupling. A convenient gauge choice, which as- 
sures a constant uniform magnetic field, is A^ = 
{0,—Bx'^/2,Bx^/2). The Dirac equation takes now the 
form: 

[i7°((9o - Jm) + ivFl^Di + ivFl'^D2\ * = (20) 
where we define: 



Di^di 



ieBx 



D2=d2 + 



ieBx 



(21) 



Inserting ^ = exp{—iEx^)^ to (PU]) leads to the following 
Hamiltonian form of the Dirac equation: 

{E + ij)<^{x) = -iVF {l^l^Di + ^°-i^D2) <^{x) (22) 



In terms of the two-component spinors <i> = 
equation (P^ reads 

{E + ii)x = -ti'F {(j'^Dk) X 
{E + h)k = ivp (a'^Dk) k 



the 



(23) 



This is the set of decoupled equations for upper and lower 
components. Thus, the Dirac Hamiltonian has the form: 



H = 



H+ 
iJ_ 



with 



H± = ^ivp 



di 



di 



^^.2. 



- id2 



eB_ 1 

2 



«^x2. 



' 182 



(24) 



£S 1 
2 



(25) 

Below we will follow the discussion given in Reffigl. It is 
useful to introduce the complex variable: z = x"^ + ix^ 
and to define the operators: 

I-- = -2| + f5 ,26) 

It can be easily found that they obey the algebra: 

[D,D*]^2eB (27) 

The angular momentum operator expressed in terms of 
the above defined complex variable reads: 



dz dz 2 



(28) 



The a^ - contribution is the so-called lattice spin and not 
the real electron spin (cf. Ref[3 and the discussion in Sec. 
IV below.) 

The Hamiltonian H± acquire a very simple form: 



H± = TVF 



D 
D* 



(29) 



This shows that the eigenvalue problem exhibits 
supersymmetrjii^. Therefore, proceeding in the standard 
way we define the following vectors: 



where ^0 obeys: 



vi'„ = iD*r^o 



D'i'o = 



(30) 

(31) 



It is now straightforward to check that $„ obey the eigen- 
value equations: 



D*D-^r. 



2neB^n 
2(nH-l)eS*„ 



(32) 



It can be easily found that the eigenvalues and the eigen- 
vectors of i7+ have the following form: 



H+Xn+ = VFV2neBxn+, Xn+ = 






H+Xn- = -UF\/2(71 + l)eBxn-, Xn- 



^2(n + l)eS 

(33) 



where n = 0,1,2,.... It should be noted that the spec- 
trum of _ff-|_ is infinitely degenerate. Namely, H^ com- 
mutes with angular momentum J, so a given eigenspace 
of 7J-|_ is spanned by the eigenvectors of J corresponding 
to arbitrary allowed values of angular momentum. To ex- 
plain this, one can observe that the equation (pijl defines 
^0 up to the arbitrary multiplicative factor depending on 
z. The general solution of the equation defining ^Pq can 
be written as a linear combination of following eigenvec- 
tors: 



*'n 



(34) 



They carry the orbital angular momentum —m and al- 
lows us to define, with the help of eq. ((30|) the towers of 
eigenvectors of DD* and D*D: 



^m ^ (£)*)"*J," 



(35) 



^Jf are eigenvectors of orbital angular momentum oper- 
ator corresponding to the eigenvalue n~m. Accordingly, 
the complete set of eigenvectors of H+ reads: 



Xn+ ^ 



V2^^^*™_i 



An — 



^ r. 



^2(n+l)eB 



(36) 



The total angular momentum carried by x™± equals 
n — r?T, =p 1/2. For H- one similarly obtain the follow- 
ing eigenvalues and eigenvectors: 



H^kI 



-Vp 



V2neB<L, k™_ 



iJ_ kI\ = VF \/2{n + l)eSC+ , k. 



■ri+ 



V2^^*™_i 

^ n 
^ n 



i/2(n+l)ei 



(37) 



The total angular momentum carried by k™^ equals 
n — 777. =p 1/2. 



Let us consider the subspace corresponding to definite 
energy, £"„ = vp \/2neB — jj, and angular momentum j — 
n — m ~ 1/2. Dirac equation takes the form: 



L 



X 







En + /i - 


-H^ 




X 


Y 




En +1-1- 


-H+ 






Y 





(38) 



Two eigenvectors corresponding to the above energy and 
angular momentum have the following form: 



62 



n 






64 



" 





^ n — 1 


*71 


L V2neB J 



(39) 



We define two further vectors forming together with 62 
and 64 a basis: 



61 



-V[/l 



*" 



2\/2neB 



63 



*" 



2V2neB 
4nei3 






(40) 



One can observe that 61, 62, 63, 64 span an invariant sub- 
space for the operator L in which L takes the form given 
by PH]) . Therefore the symmetry algebra of Dirac equa- 
tion restricted to the subspace of definite energy and an- 
gular momentum is: 



S = {gL{2, C) ® gL{2, C))+ A{4, C) 



(41) 



We have found that the symmetry algebra coincides 
with the one corresponding to the free case. However, 
both cases differ by the choice of commuting variables 
defining the relevant subspace. In the free case these 
were the momentum components while here we consider 
the subspaces of given energy and angular momentum. 



IV. DISCUSSION AND PERSPECTIVES 



by different choices of commuting operators (the momen- 
tum components in the free case and energy and angular 
momentum in the case of uniform magnetic field). 

The main difference between the standard 3+1- 
dimensional Dirac equation and the one considered here 
is that in the former the 7 matrices form the (unique up 
to equivalence) irreducible representation of Clifford alge- 
bra. On the contrary, in the 2+1-dimensional case there 
are two, both two-dimensional, inequivalent irreducible 
representations. When writing out the Dirac equation 
we use both of them forming the direct sum which pro- 
vides four-dimensional reducible representation. As a re- 
sult, not all symmetry operators can be constructed out 
of 7 matrices (and their products) entering the relevant 
Dirac equation; one has to use, additionally, the 73 ma- 
trix which is absent from the relevant Dirac equation. 

Let us note that the matrix structure which determines 
the rotational symmetry has all properties of spin alge- 
bra, in spite of the fact that the internal degrees of free- 
dom come from the existence of two sublattices rather 
than from electron spin degrees of freedom which are 
neglected in the considered approximation. In fact, it 
has been argued^ that the pseudospin, arising from the 
degeneracy introduced by the honeycomb lattice's two 
atomic sites per unit cell, has properties of real angular 
momentum. The form of the solutions of Dirac equation 
in magnetic field which we have used to find the symme- 
try algebra supports this point of view: we are dealing 
with the solution describing spin-i particles in magnetic 
field. The problem of "lattice spin" seems to be worth 
further considerations. In high energy physics the emer- 
gence of spin from other degrees of freedom is known 
phenomenon called "spin from isospin"— 1^. For exam- 
ple, one can consider isospin gauge theory with bosonic 
matter of isospin i (as well as neutral Higgs boson). It 
possesses the monopole solution. Quantizing the the- 
ory in one- monopole sector one finds that, due to the 
very presence of the monopole the SU(2)rotationsxSU(2) 
isospin symmetry is broken down to its diagonal sub- 
group which is the symmetry of monopole configuration. 
The relevant generators are the sums of ordinary rota- 
tions and isorotations. As a result, with bosons of half- 
integer isospin one obtains half-integer spin, in spite of 
the fact, that there are no elementary fermions involved. 
It would be interesting to interpret the results of Ref|^ 
in a similar spirit. 



We have analyzed the "on-shell" symmetries of 2+1 
dimensional Dirac equation for the free particle as well 
as the one interacting with uniform magnetic field. In 
both cases the symmetry algebra appears to be the same 
(cf. eq. ([41]) ). However, it refers to the subspaces defined 
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